ARL  67-0187 
SEPTEMBER  1967 


Aerospace 


Laboratories 


A  CRITICAL  REVIEW  OF 
HETEROGENEOUS  MIXING  PROBLEMS 


ANTONIO  FERRI 
NEW  YORK  UNIVERSITY 
BRONX,  NEW  YORK 


Conlracl  No.  AF  33(615)-2215 
Project  No.  7064 


This  document  has  been  approved  for  public  release  and  sale; 
its  distribution  is  unlimited. 


DISCLAIMER 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


I 


»  AsfiouHcro 
atiific  ;im 


s'.M  ^ 

[J 


NOTICKS 


ode  mitwitirr  rrTil 


<iM(  YVherf.’Gtti  eminent  drawings,  specifications,  or  other  data  are  iim  iI  for  any  purpose  otliei  than  in 
'  connection  wit  i  a  definitely  related  Government  procurement  operation,  the  I’nited  States  Can  eminent 
thereby  incurs  no  responsibility  nor  any  obligation  wliutsoevoi ,  and  the  fact  that  the  Government  may 
have  f  irnuilatc  1,  furnished.  or  in  any  way  supplied  tin4  said  drawings,  specifications,  or  other  data,  is 
not  to  be  rega  (led  by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any  other 
person  or  corporation,  or  conveying  any  rights  or  piinussion  to  manufacture,  use,  or  sell  iny  patented 
invention  that  may  in  any  way  be  related  thereto. 


Agencies  of  the  Department  of  Defense,  qualified  contractors  and  other 
government  agencies  may  obtain  copies  from  the 

Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  22314 

This  document  has  been  released  to  the 


CLEARINGHOUSE 

U.S.  Department  of  Commerce 

Springfield,  Virginia  22  151 


for  sale  to  the  public. 


Copies  of  ARL  Technical  Documentary  Reports  should  not  be  returned  to  Aeiospace  Research 
Laboi atones  unless  return  is  required  by  security  considerations,  contractual  obligations  or  notices  on 
a  specified  document. 


400  -  November  1967  -  C045S  -  13-270 


IB 


ARL  67-0187 


A  CRITICAL  REVIEW  OF 
HETEROGENEOUS  MIXING  PROBIEMS 


ANTONIO  FERRI 

NEW  YORK  UNIVERSITY 
BRONX,  NEW  YORK 


SEPTEMBER  1967 


Contract  AF  33(615)  2215 
Project  7064 


This  document  has  been  approved  for  public  release 
and  sale;  its  distribution  is  unlimited. 


AEROSPACE  RESEARCH  LABORATORIES 
OFFICE  OF  AEROSPACE  RESEARCH 
UNITED  STATES  AIR  FORCE 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO 


FG2.2W0RD 


This  interim  report  was  prepared  by  Dr.  Antonio  Ferri,  Director, 
N.Y.U.  Aerospace  Laboratories,  and  presents  research  carried  out  under 
Contract  No.  AF33(615)2215,  "Boundary  Layer  Characteristics  for  Hypersonic 
Flow  in  the  Presence  of  Mass  Addition",  Project  No.  7064. 


ii 


ABSTRACT 


A  critical  review  of  mixing  problems  is  presented.  The  limitations 
of  boundary  layer  types  of  analyses  are  discussed.  A  method  of  analysis 
which  does  not  use  boundary  layer  approximation  is  outlined.  Difficulties 
of  correlating  experimental  data  attributable  to  not-well-def ined  boundary 
conditions  are  described.  New  phenomena  related  to  mixing  with  pressure 
gradient  and  chemical  reaction  are  presented. 
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1.  INTRODUCTION 


Thc  problc"  °f  ”u,n8  "»  ««—  *■  „,le,  Cla.sl<:a 
aerodynamics  has  Investigated  tot  many  year,  the  problem  of  mixing  o£  t„„ 

K  .. t reams  However,  the  quantitative  understanding  of  such  phy. 
MC  Phenomenon  Is  sell,  Incomplete  even  for  „x,„g  uhcre 

properties  are  known,  and  the  progress  made  both  by  theoretical 
»n<l  experiments  1  Investigation  Is  st,U  not  even  satisfactory. 

f'“  “  an•1lVtlCa,  °f  »'"•  ‘"e  difficulty  otlses  hecouse  of 

'"C  “"PU'“U°"  ‘»  equations.  I„  addition,  the  ,ach  of 

knowlcdpe  of  transport  properties  Increases  the  difficulties  of  the  prob¬ 
lem  for  turbulent  mixing.  Classically,  mixing  problems  are  analysed  by 
introducing  al,  simplifications  used  In  boundary  layer  theory.  u„f0r- 
tunately.  the  Introduction  of  such  simplifications  mixing  problems  are 
net  necessarily  lust, fled  and  acceptable  for  the  determination  of  Important 
Pliymca,  properties  of  the  phenomena.  Several  questions  can  be  raised  with 
respect  to  the  use  of  the  boundary  layer  equations  mixing  problems  *, 
introduction  o,  the  boundary  layer  equations  In  the  analysis  of  flc  near 

*  P""‘tS  “*«■  momentum  equation.  ,t  can  be  shown 

. .  elm" 1 1  flea t Ion  I.  not  Important  In  the  determination  of  the  var- 

mt.on  of  momentum  along  streamlines,  and  therefore  of  the  shin  friction 

“  °"C  tl,C  '•  *  tlctcrmlncd  It,  boundary  layer 

analysis.  However,  In  mixing  problem,  the  most  Interesting  Infect, on 

requ  from  the  analysis  Is  thc  rate  of  diffusion  o,  a  fluid  of  one  screen 

‘mo  the  fluid  of  the  other,  and  such  motion  Is  controlled  by  thc  second 
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momentum  equation.  As  a  consequence,  the  use  of  boundary  layer  equations 
can  be  Justified  only  when  the  characteristics  of  the  two  streams 
(velocity,  density,  composition-temperature )  are  not  too  different.  If 
the  two  streams  have  physical  properties  which  are  widely  different,  then 
the  scale  for  the  determination  of  the  order  of  magnitude,  which  is  the 
thickness  of  the  layer  affected  by  the  flow,  becomes  extremely  large. 

Typical  Indication  of  the  Inadequacy  of  the  approximation  of  the  boundary 
layer  theory  is  given  by  the  results  of  the  analysis  of  a  jet  discharging 
in  a  quiescent  atmosphere  made  with  this  approximation.  In  this  case,  in 
the  external  flow  the  velocity  component  normal  to  the  axis  of  the  Jets 
are  larger  than  the  tangential  component  in  the  proximity  of  the  origin  of 
mixing,  therefore,  the  description  of  the  motion  of  the  external  flow  can¬ 
not  be  represented  by  the  simplified  equations.  The  normal  component  of  the 
velocity  is  directly  related  to  the  mass  of  the  flow  entrained  by  the  jet, 
and  therefore  is  a  quantity  that  must  be  determined  with  satisfactory 
precision. 

When  the  approximation  of  the  boundary  layer  theory  cannot  be 
Justified  for  the  analysis  of  mixing  problems,  then,  as  a  consequence, 
the  boundary  conditions  for  a  mixing  problem  cannot  be  specified  in  a 
simple  form  at  the  station  where  mixing  begins  if  the  two  streams  are 
subsonic.  In  this  case,  the  mixing  process  induces  velocity  components 
which  are  significant  upstream  of  the  initial  mixing  region. 

For  the  same  reason  it  is  very  difficult  to  perform  significant 
exncriments  of  mixing  when  the  two  streams  are  quite  different 
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In  chla  case,  the  characteristics  of  the  experimental  apparatus 
upstream  of  the  mixing  region  have  important  effects  on  the  mixing 
phenomena  near  the  origin  of  mixing.  This  point  is  often  neglected  in 
presenting  experimental  reaults  and  in  comparing  experimental  results. 

When  the  two  streams  are  subsonic  even  for  the  case  of  the  mixing  of  two 
low  velocity  streams  of  the  same  fluid,  the  initial  conditions  cannot 
be  defined  by  giving  the  average  properties  of  the  two  streams  as,  for  cxamp 
the  average  velocities  measured  in  terms  of  stream  tube  area  and  mass 
flow,  but  are  defined  only  if  the  complete  velocity  distribution  and 
pressure  distribution  at  the  station  at  the  beginning  of  the  mixing 
are  given.  Such  distributions  can  change  significantly  from  experiment 
to  experiment  because  they  are  affected  by  the  experimental  apparatus 
upstream  of  the  mixing  region;  therefore,  presentation  of  experimental 
data  on  the  basis  of  simple  parameters  as  average  velocity  ratio  or 
average  velocity  difference  of  the  two  jets  j.s  not  Justified  because 
such  terms  do  not  define  cither  physical  properties  or  the  experiment. 

2.  ANALYSIS  OF  COMPLEX  VISCOUS  FLOW  PROBLEMS 
When  the  properties  of  the  two  mixing  streams  arc  quite  different 
the  solution  of  more  complex  sets  of  equations  is  required  where  pressure 
gradients  in  two  perpendicular  directions  are  considered  to  be  of  the 
same  order  and  where  viscous  effects  arc  considered.  Such  type  of 
analysis  is  of  great  impoitsnce  in  many  problems  of  practical  interest 
besides  the  chemical  problem  of  coaxial  mixing  of  nonreacting  streams. 
Typical  examples  of  phenomena  where  such  type  of  analysis  is  required 


are:  (a)  The  mixing  of  two  reacting  streams  where  the  energy  released 

by  the  reaction  produces  important  pressure  perturbation  that  cannot 
be  neglected;  (b)  High  speed,  low  density  flow  around  bodies;  in  a 
wide  Reynolds  number  range,  for  many  practical  shapes  the  flow  still 
can  be  analyzed  with  the  assumption  of  continuous  medium,  but  gradients 
of  velocity  arc  important  in  the  flow  fields  and  therefore  important 
viscous  effects  are  present  in  all  flow  fields;  (c)  Boundary  layer  with 
large  rates  of  ablation  where  the  blowing  velocity  is  large;  (d)  Boundary 
layer  with  chemical  reactions  where  the  chemical  enthalpy  of  reaction 
is  of  the  same  order  of  the  total  enthalpy  of  the  flow;  (e)  Hypersonic 
boundary  layers  with  large  pressure  gradients  along  the  streamline. 

For  such  problems  the  variation  of  pressure  in  the  boundary  layer 
thickness  is  of  the  same  order  of  the  variation  between  two  distant 
points  along  a  given  streamline  because  M  is  large.  In  this  case,  the 
static  pressure  varies  strongly  in  boundary  layer  in  normal  direction 
to  the  wall  and  cannot  be  neglected.  The  effect  of  such  pressure 
gradient  affects  strongly  the  velocity  profile  and  therefore  all  the 
viscous  effects.  An  example  of  such  pressure  variation  is  shown  in 
Fig.  1.  The  upper  part  of  the  figure  shows  th  •  experimental  apparatus, 
while  the  lower  part  shows  the  pressure  varlat  on  and  velocity  distribution 
in  the  boundary  layer  at  a  station  of  the  body.  The  pressure  at  that 
station  is  27  tiroes  the  initial  static  pressure.  The  flow  Is  initially 
at  M  *  6,  and  the  body  Is  axially  symmetric.  The  boundary  layer  thickness 
is  of  the  order  of  0.35".  The  experiments  have  been  performed  at  the 
Laboratory  of  New  York  University  by  Mr.  Hoydysh  (Doctoral  candidate). 


The  solution  o i  viscous  problems  where  more  complex  equations  are 
required  can  be  approached  along  two  different  lines.  The  classical 
line  Is  the  perturbation  method  where  higher  order  terms  are  obtained 
as  second  approximation.  The  calculation  starts  with  the  boundary 
layer  approximation  which  is  assumed  as  first  order  solution  In  the 
determination  of  the  flow.  This  method  permits  us  to  obtain,  In  some 
Instances, a  second  order  correction;  however,  such  an  approach  Is  use¬ 
ful  only  If  the  first  order  solution  gives  a  sufficiently  good  approximation 
of  the  flow  field,  then  the  second  approximation  can  be  used  to  improve 
the  accuracy  of  the  results.  For  many  of  the  problems  described  such 
an  approach  is  not  very  promising. 

A  second  approach  suggested  by  the  author  and  Dr.  Morcttl1  consists 
in  simplifying  the  equation  of  motion  on  the  basis  of  the  following 
physical  principle.  The  viscous  effects  arc  produced  by  gradients. 

In  many  physical  problems  such  as  those  previously  described,  the 
variation  of  physical  properties  along  streamlines  is  much  smaller  than 
in  a  direction  normal  to  the  streamlines,  then  the  cquat ions  that  govern 
the  motion  of  the  fluid  can  be  written  with  good  approximation  in  the 
following  form: 

In  the  two-dimensional  case,  the  momentum  equations  for  the  above 
formulation  arc: 


Oqq9  +  Ps 

“  A, 

(1) 

pq309  +  Pn 

-  o. 

U) 

where  p  is  density,  p  is  pressure,  q  and  6  the  modulus  and  the  direction 
of  the  velocity  relative  to  engine  axis.  The  quantity  A  is  given  by 
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A '  HV 


where  ^  is  the  viscosity. 


The  subsequent  solution  of  the  compatibility  equations  requires 

that  A  be  constant  or  of  small  relative  variation  during  an  integration 

step.  Then  the  initial  computation  can  be  made  approximately  using  the 

initially  known  value  A°,  followed  by  iteration  to  obtain  a  more  accurate 
mean  value  A. 

The  continuity  equation  reads 

qps  +  pqa  +  pq0n  «=  0.  <4) 

and  the  energy  equation  can  be  stated  (assuming  for  simplicity  that  the 
Lewis  and  Prandtl  numbers  are  equal  to  unity) 


where  H  is  the  total  enthalpy 


h  being  the  static  enthalpy. 


H  -  — (uH  )  . 
s  pq^n'n* 


H  "  h  +  f-  , 


The  diffusion  along  the  streamline  has  been  neglected  in  comparison 


with  lateral  diffusion. 


In  order  to  complete  the  description  of  the  system  we  express  the 
enthalpy  by  means  of  the  equation 


h  "  )  ha  +  T  >  Q  — i 

"  <—  l  is  s  L  i  dt 

l°l  i-1 


where  Qt  is  the  mass  fraction  of  species  i.  From  the  equation  of  state 


we  obtain  further 


PC3  PTS  V  QiS 

• " — r  -  ^  L  -$r  ■  °* 
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R  being  the  gas  constant  and  the  molecular  weight  of  species  i. 

The  species  conservation  equations  must  take  into  account  the  mixing 
of  species  between  streamlines  and  the  change  in  mass  fraction  of  the 
species  due  to  chemical  reaction  along  the  streamline.  Thus, 


i  a  Vi 

a.  ■  —  t-  pD  t —  +  - 

is  pq  on  on  pq 


(9) 


where,  again  for  simplicity,  a  single  binary  diffusion  coefficient  D 
has  been  assumed.  For  unit  Lewis  number  we  can  substitute 

pD  -  p  (10) 

The  subsequent  algebraic  manipulations  to  put  the  equations  into 
suitable  characteristic  form  arc  tedious,  so  details  should  be  sought 
in  Ref.  6.  It  is  shown  that  the  third  equation  to  be  used  with  (1) 
and  (2)  can  be  put  in  the  form 


(q  +  -)  q  +  -  p  +  f6  -  e, 

v  qj  Ms  p  rs  n 


(11) 


where 


dh, 


F  "  T  2.  Q1  It  ' 


E  -  —  (pH  )  -  j  h  a  +  ,  r~  . 

pq  'n  n  U  11  p  W 


a. 


1  1  * 

It  is  seen  that  the  viscosity  docs  not  appear  in  (l),  (2)  and  (11) 
except  n  the  right-hand  side  and  so  does  not  directly  affect  the  streamline 
and  characteristic  directlonr  and  the  location  of  the  net  Intersections. 

The  quantities  related  to  transport  properties  arc  represented  by 
the  terms  in  the  right-hand  side  of  Eqs.  (1),  (2)  and  (11).  Such 
quantities  arc  functions  of  the  coordinates,  and  of  the  gradients  of 
local  flow  properties  and  therefore  can  be  determined  only  by  the 
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solution  of  the  differential  equations.  However,  in  practical  numerical 

solutions,  the  values  of  such  terms  can  be  estimated  with  sufficient 

precision  from  the  values  of  the  same  quantities  at  neighboring  points. 

On  the  basis  of  such  an  estimate  a  numerical  value  can  be  determined 

for  the  local  value  of  such  quantities.  Then  the  left-hand  side  of  the 

equations  can  be  solved  by  using  known  methods  of  analysis  for  inviscid 

flows.  When  the  flow  is  determined  in  first  approximation,  a  second 

approximation  is  feasible  where  these  terms  in  the  right-hand  side  of 

Eqs.  (1),  (2)  and  (11)  are  calculated  on  the  basis  of  the  local  flow 

properties  determined  from  the  first  approximation.  Such  an  approach 

is  very  useful  when  the  "inviscid"  part  of  the  flow  contributes  appreciably 

to  the  local  variation  of  flow  properties,  while  the  effect  of  transport 

phenomena  is  smaller  because  it  is  more  gradual.  In  this  case  boundary 

layer  type  of  equations  cannot  give  a  satisfactory  first  order  approximation. 

I  2 

This  approach  has  been  used  for  the  case  of  supersonic  t  lows  ’  .  In 
this  case,  the  inviscid  part  is  solved  with  the  method  of  characteristics 
by  means  of  a  difference  method,  by  determining  as  a  first  step  the 
effect  of  terms  containing  the  transport  properties  at  the  two  initial 
points,  l  and  2  of  Fig.  2,  where  the  properties  arc  known.  Such  terms 
arc  reevaluated  in  second  approximation  at  point  3,  and  average  values 
between  1,  3  and  2  and  3  arc  used.  Dr.  Moretti  has  applied  such  an 
approach  to  problems  of  supersonic  mixing  with  chemical  reaction. 

Results  of  calculations  for  such  problems  are  shown  in  Figs.  3  and  U. 

The  problem  analyzed  is  mixing  with  combustion  of  two  coaxial  axially 
symmetric  streams  of  hydrogen  and  air.  The  two  jets  are  initially 
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parallel  and  at  the  same  static  pressure;  however,  the  mixing  and 
combustion  pioduce  waves  that  propagate  in  the  flow  field.  In  Fig.  3 
the  isobars  are  shown.  The  combustion  produces  pressure  waves  that 
propagate  in  the  flow  outside  of  the  mixing  region.  In  Fig.  4,  the 
isotherms  for  the  same  case  are  shown.  The  pressure  field  outside 
of  the  mixing  region  is  obtained  as  a  result  of  the  malysis.  If  the 
same  problem  is  analyzed  by  means  of  classical  mixing  type  of  analyses 
based  on  boundary  layer  equations,  then  the  pressure  field  outside  of 
the  mixing  region  must  be  assumed  as  given  and  cannot  be  determined 
analytically. 

Another  application  of  the  approach  described  is  related  to  the 
determination  of  the  low  density  supersonic  flow  around  bodies  as  shown 
in  Fig.  5.  The  Reynolds  number  is  small,  therefore  the  displacement  of 
the  boundary  layer  affects  the  pressure  downstream  of  the  nose.  In 
addition,  because  of  the  curvature  of  the  shock,  the  velocity  gradients 
in  the  flow  field  are  finite  and  the  shock  thickness  cannot  be  neglected, 
therefore  the  viscosity  is  important  in  all  flow  fields.  Such  types 
of  flow  are  usually  analyzed  by  neglecting  the  viscosity  between  the 
shock  and  the  body  surface,  and  by  determining  the  viscous  effects 
at  the  wall  with  boundary  layer  analysis  in  the  presence  of  external 
vorticlty.  Such  an  approach  requires  some  discontinuity  between  the 
two  flow  regions  independently  of  how  the  analysis  is  perform.,  d.  The 
discontinuity  in  assumptions  of  the  two  flow  fields  is  not  Justified 
due  to  the  presence  of  curvature  of  the  shock.  The  transport  properties 
can  be  considered  in  all  flow  fields  if  the  flow  is  analyzed  with  the 


9 


approach  described  before.  In  a  small  layer  near  the  wall,  where  the 
curvature  of  the  streamline  is  small,  boundary  layer  types  of  analyses 
are  used.  The  flow  at  the  outer  edge  is  supersonic.  Such  a  layer  is 
much  smaller  than  the  boundary  layer,  therefore,  boundary  layer  type  of 
equations  are  sufficiently  accurate.  The  matching  with  the  external 
flow  takes  into  account  quantities  and  derivatives  both  for  velocity 
and  enthalpy.  This  can  be  done  because  transport  properties  also 
are  considered  in  the  equations  that  define  the  external  flow.  The 
external  flow  is  analyzed  by  means  of  the  method  of  characteristics 
in  presence  of  transport  effects  as  described  previously.  The  con¬ 
ditions  at  the  shock  take  into  account  the  effect  of  shock  thickness 
as  discussed  in  the  doctoral  thesis  of  R.  R.  Chow'*,  (this  work  has 
been  performed  by  Dr.  E.  Krause  at  NYU  as  his  doctoral  thesis).  The 
basic  equations  for  the  three  regions  are  indicated  in  Figs.  5-7. 

Such  type  of  analysis  can  be  applied  in  the  same  way  to  the  analysis 
of  boundary  layer  with  large  pressure  gradients,  when  the  transport 
properties  are  known. 

The  flow  properties  near  the  wall  are  determined  with  boundary 
layer  equations  as  shown  in  Fig.  6,  while  the  outer  flow  which  is 
supersonic  is  analyzed  with  the  equations  of  Fig.  7. 

3.  TURBULENT  MIXING 

In  addition  to  the  difficulties  of  an  analytic  nature,  additional 
difficulties  are  present  in  the  case  of  turbulent  mixing  due  to  lack 
of  knowledge  of  the  physical  laws  that  regulate  the  transport  properties 
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in  turbulent  flows.  The  transport  properties  for  turbulent  flow  are 
affected  by  many  different  characteristics  of  the  flow;  at  present, 
even  the  most  Important  parameters  which  affect  such  l  msport  properties 
are  not  completely  determined.  Because  of  the  complexity  of  the 
phenomenon,  it  is  traditional  to  attempt  to  relate  these  transport 
properties  to  a  few  of  the  global  flow  properties  and  to  determine 
their  values  experimentally  as  a  function  of  such  parameters.  However, 
the  approximate  dependence  of  such  properties  on  gross  values  of  the 
flow  parameter  is  not  well  understood.  The  reasons  are  several.  Such 
dependence  is  obtained  on  the  basis  of  experiments.  The  experimental 
work  available  is  limited,  and  many  of  the  experiments  are  not 
satisfactorily  defined  for  such  a  case.  In  many  of  the  experiments 
even  the  statement  that  the  mixing  is  turbulent  must  be  carefully 
scrutinized.  In  boundary  layer  investigations  the  turbulent  boundary 
layer  has  a  region  of  laminar  motion  near  the  wall  while  at  some 
distance  from  the  wall  random  turbulent  fluctuations  exist  and  farther 
on  they  disappear.  The  persistence  of  turbulent  motion  is  related  to 
a  Reynolds  number  for  which  a  characteristic  velocity,  density,  and 
length  that  defines  the  experiment  must  be  defined. 

In  a  mixing  of  two  coaxial  Jets,  two  characteristic  velocities 
exist,  two  different  densities  must  be  considered,  and  a  characteristic 
length  is  difficult  to  define.  For  some  flow  conditions  the  random 
motion  can  persist  in  one  of  the  two  streams  while  the  motion  in  the 
other  stream  can  be  essentially  of  laminar  type.  It  is  therefore 
probable  that  the  transition  for  a  purely  laminar  to  a  completely 
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turbulent  nixing  Is  very  gradual,  and  the  transit  tonal  conditions 
could  cover  a  wide  range  of  combination  of  Reynold*  numbers  of  the 
two  *  t  re  aaa . 

In  thla  case  the  detc rialnat Ion  of  preaencc  of  turbulence,  for 
example,  by  Dean*  of  a  hot  wire  anemometer,  In  the  high  velocity  stream 
does  not  neccaaarlly  define  a  purely  turbulent  nixing,  but  could 
Indicate  a  transitional  region.  For  example.  It  la  probable  that  the 
type  of  phenomena  obtained  when  both  streams  have  large  velocity  and 
density,  and  therefore  high  Reynolds  numbers,  are  different  from  the 
profiles  when  the  streams  have  a  low  density  and  the  satoc  velocity 
as  In  the  first  case,  even  If  one  stream  Is  already  turbulent  while 
the  other  Is  Initially  laminar.  The  problem  Is  more  difficult  when 
one  of  the  streams  has  very  low  velocity  and  therefore  very  low  Reynolds 
number. 

A  second  problem  that  makes  '  t  difficult  to  determine  the  dependence 
of  the  transport  properties  on  gross  flow  properties  la  related  to  the 
difficulty  of  accurately  defining  the  initial  conditions  and  boundary 
conditions.  Usually  the  properties  of  the  two  streams  are  defined 
in  terms  of  average  properties  at  the  initial  section*.  Such  properties 
arc:  chemical  species  of  the  flows,  enthalpies,  and  velocities.  Often, 
average  values  of  such  quantities  nrc  assumed  to  be  representative,  which 
corresponds  to  step  function  distributions  for  the  three  quantities. 

In  addition. the  static  pressure  external  to  the  mixing  Is  assumed  to 
be  practically  constant.  However,  the  actual  conditions  arc  much  more 
complex,  especially  for  subsonic  mixing.  The  mixing,  especially  when 
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the  velocltrt  of  the  two  jots  are  very  different,  tends  to  produce 
larj'i’  longitudinal  and  radial  pressure  gradients  and  therefore  large 
voloclly  gradients  are  present  In  the  region  near  the  origin  of  nixing. 
Such  eflects  travel  upstream  and  Induce  a  flow  field  lar  flora  mil  forra, 
therefore  the  flow  and  properties  at  the  beginning  of  mixing  ale  far 
from  the  average  values  assumed,  therefore  the  gradients,  which  are 
related  to  the  transport  phenomena,  In  this  region  cannot  be  defined 
only  by  stop! e  average  parameters.  In  addition,  the  presence  of 
boundary  1  avers  mod  1 f y  substantially  the  average  distribution  both  of 
velocity  and  enthalpy  The  heat  transfer  across  the  wall  that  divides 
the  coaxial  jet#  can  have  serious  effects  on  the  flow  in  the  near  re¬ 
gion  ol  mixing  In  experiments  where  the  enthalpies  of  the  two  streams 
are  different  These  effects  are  reduced  when  the  external  stream  Is 
supersonic  and  the  axial  stream  has  a  cylindrical  entrance  region.  In 
this  case  the  l 1 ow  at  the  beginning  of  the  mixing  can  be  defined  fairly 
accurately,  Independently  of  the  mixing  effects.  In  subsonic  flow, 
an  additional  problem  must  he  considered  which  is  related  to  the  Inter¬ 
ference  of  the  boundary  of  the  external  jet  an  Indicated  In  Dr, 

l. 

Abramovich  In  his  paper 

In  subsonic  flow  the  turbulence  generated  at  the  bout dary  of  the 
external  jet  can  propagate  and  affect  the  mixing  region.  Again,  such 
an  effect  can  be  eliminated  in  supcrFonlc  1 1 ow  if  the  mixing  takes 
place  upstream  of  the  Kach  cone  from  the  boundary.  (Tills  was  the  case 
In  the  experiments  performed  under  my  supervision,  contrary  to  what 
has  been  stated  In  Ref.  A. )  The  presence  of  Induced  flow  In  subsonic 


streams  Is  strongly  dependent  on  the  boundary  conditions,  and  therefore, 
for  subsonic  flow  It  Is  more  difficult  to  obtain  Information  on  the 


basis  of  simple  parameters.  For  example,  when  the  velocity  of  the 
external  stream  1s  ouch  smaller  than  the  velocit>  of  the  other  Jet,  It 
seems  arbitrary  o  try  to  define  the  properties  of  mixing  on  the  basis 
of  quantities  related  to  two  average  velocities,  while  the  velocities 
In  the  lczscdlatc  region  outside  of  the  mixing  arc  quite  different  than 
such  values. 

The  quantity  that  Is  less  affected  by  the  flow  field  distribution 
and  by  the  boundary  conditions  Is  the  cor- *nt rat  Ion  of  species  which 
Is  given  at  the  beginning  of  mixing  always  as  a  step  function,  there¬ 
fore  in  heterogeneous  mixing  this  quantity  Is  more  uniformly  representative 
of  the  initial  conditions. 

4,  PRESSURE  GRADIENT  EFFECTS 

From  an  experimental  point  of  view,  axially  symmetric  Jets  are 
usually  preferred  to  "two-dimensional"  type  of  Jets  because  of  the 
difficulty  of  producing  really  two-d Imcns  ional  Jets  with  finite  side 
boundaries.  However  in  axially  synxaetric  supersonic  Jets  problems  of 
other  nature  arc  present  that  often  arc  neglected.  Adverse  pressure 
gradients  in  a  supersonic  stream  arc  transmitted  by  waves.  In  axially 
syrraotrlc  flow,  waves  become  stronger  at  the  axis  of  the  Jet.  Such 
focusing  of  waves  produces  much  larger  pressure  gradients  at  the  axis 
than  at  some  distance  of  the  axis;  such  radial  pressure  variation 
strongly  affects  the  velocity  gradients  because  it  changes  the  radial 
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velocity  profile,  and  therefore  affects  the  mixing.  Small  gradients 
usually  are  always  present  In  wind  tunnels,  or  in  wakes,  even  at  some 
distance  of  the  body  especially  If  the  body  travels  In  supersonic 
stream.  Such  small  gradients  can  strongly  affect  the  mixing  process 
and  cannot  be  neglected  as  has  been  usually  done.  In  order  to  Indicate 
the  Importance  of  such  affects  an  experiment  has  been  performed,  at 
my  request,  by  Dr.  Zakkay  at  New  York  University  Laboratory.  A  wake 
has  been  produced  experimentally  at  M  “  3.89,  at  high  Reynolds  numbers. 

The  Mach  number  along  the  axis  In  a  region  of  the  wake  has  been  determined 
for  the  case  of  external  pressure  practically  uniform  oc  the  wall  of  the 
tunnel.  In  a  second  test,  a  strip  of  the  thickness  of  the  order  of 
1/20  inn  (scotch  tape)  has  been  glued  at  the  wail  of  the  tunnel.  The 
measurements  have  been  repeated.  The  difference  In  results  between 
the  two  experiments  Is  shown  In  Fig.  8.  The  pressure  disturbance 
produced  by  the  tape  at  some  distance  from  the  axis  Is  smaLl  and  ol 
the  order  of  precision  of  experiments.  However,  the  effect  at  the  axis 
Ib  significant.  The  effects  of  pressure  gradients  can  be  very  large 
In  mixing  phenomena,  and  can  modify  substantially  the  flow  field 
properties  of  the  flow. 

5.  MIXING  WITH  ADVERSE  PRESSURE  GRADIENTS 
The  presence  of  longitudinal  pressure  gradients  is  of  extreme 
Interest  because  It  can  produce  flow  fields  which  are  quite  different 
from  those  corresponding  to  zero  pressure  gradients.  Consider  the 
mixing  of  two  heterogeneous  coaxial  streams.  Assume  that  an  adverse 
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pressure  gradient  la  Imposed  In  the  outside  flow  and  Is  extended  to  a 
length  such  that  the  total  pressure  rise  can  produce  reverse  flow  In 
the  stream  at  lower  Mach  number  If  mixing  Is  neglected.  The  mixing 
can  accelerate  the  lower  Mach  number  stream,  and  therefore , the  ability 
of  such  a  stream  to  traverse  the  pressure  gradient  depends  on  the 
amount  of  shear  force  produced  by  the  mixing  process.  In  Fig.  9  the 
result  of  such  analysis  Is  shown  for  two  assumed  values  of  eddy  viscosity, 
one  having  twice  the  value  of  the  other. 

For  the  larger  value  of  the  eddy  viscosity  the  flow  goes  through 
the  pressure  gradient  while  for  the  smaller  value  reverse  flow  occurs. 

Figure  10  gives  the  Mach  number  profiles  at  several  stations  for 
the  two  cases. 

Such  an  investigation  can  be  used  for  the  experimental  determination 
of  average  values  of  turbulent  eddy  viscosity. 

The  presence  of  adverse  pressure  gradients  in  the  presence  of 
mixing  Is  the  determining  factor  for  the  definition  of  the  flow  field 
In  the  region  of  the  base  flow  behind  a  body.  The  flow  at  the  base  of 
a  body  inrocrscd  In  an  external  stream  Is  the  limiting  caoc  of  coaxial 
mixing,  and  corresponds  to  the  case  In  which  the  momentum  of  the  central 
jet  is  zero.  If  the  Impulse  of  the  central  Jet  increases  from  zero  to 
a  finite  value  the  flow  field  changes  gradually  from  the  flow  corresponding 
to  a  base  flow  In  the  mixing  region  to  the  flow  corresponding  to  the 
mixing  of  two  parallel  streams,  both  having  large  velocity. 

The  flow  corresponding  to  this  transition  Is  very  complex.  An 
analysis  of  such  flow  indicates  the  existence  of  flow  fields  that  are 


16 


riot  well  known  at  present.  A  qualitative  desc  rlpt  lcn  obtained  from 
an  experimental  Investigation  performed  at  1«YL'  by  Dt.  Zakkay  la  presented 
In  Fig.  12  and  following  (Ref.  5).  The  experimental  apparatus  is 
shown  In  Fig.  11.  The  external  flow  is  supersonic  at  M  -  3.89.  The 
central  body  Is  a  streamline  of  the  nozzle  flow  for  uniform  flow  at 
the  beginning  of  the  test  section.  The  external  gas  Is  air,  while 
different  constituents  have  be*  n  used  in  the  central  flow.  In  Fig. 

12,  four  steps  of  the  transition  f roa  base  flow  to  a  nixing  type  of 
flow  are  represented.  In  the  upper  figure,  no  flow  Is  Injected  at  the 
bases.  The  flow  recirculates  In  the  cavity.  In  Fig.  12b  a  snail 
anount  of  flow  Is  Injected  uniformly  at  the  base  through  a  porous 
surface  placed  at  the  end  of  the  Jet.  Because  of  the  radial  and 
longitudinal  pressure  gradients  a  vortical  region  which  has  a  ring 
shape  la  produced  near  the  axis  and  the  rear  stagnation  point  raoves 
from  the  center  of  the  base  somewhat  downstream  of  tin  base.  If 
the  mass  flow  Injected  Increases,  the  recirculating  region  gradually 
moves  downstream  until  the  stagnation  pressure  of  the  injected  gas 
1s  sufficiently  high  to  overcome  the  adverse  pressure  gradient  produced 
by  the  external  flow.  Then  the  recirculating  region  moves  away  and 
disappears.  The  pressure  at  the  base  of  the  central  body  increases 
with  the  mass  flow  Injected  until  a  point  where  the  recirculating 
region  moves  downstream  and  disappears  and  then  remains  roughly 
constant.  In  Fig.  13  the  pressure  is  given  as  a  function  of  the  non- 
dimensional  izod  coefficient  for  different  gases. 

The  effect  of  pressure  gradients  becomes  more  critical  when 


17 


•  : 


chemical  reactions  take  place,  because  In  this  case  the  release  of 
energy  due  to  chemical  reaction  increases  the  speed  of  sound,  therefore 
decreases  the  Mach  number  making  easier  the  production  of  reverse  flow. 
It  is  interesting  to  observe  that  while  In  the  case  of  mixing  without 
reaction  the  Initial  point  for  reverse  flow  is  the  center  of  the  Jet, 
as  ohown  in  Fig.  14.,  in  the  case  of  chemical  reaction  the  Initial 
region  of  reverse  flow  can  be  in  the  combustion  region  that  Is  at  some 
distance  from  the  axis.  The  re  fore,  the  effect  of  pressure  gradient  for 
mixing  with  chemical  reaction  can  be  different  than  for  frozen  flows. 
Consider,  for  example,  the  mixing  problem  shown  in  Fig.  14.  A  Jet  of 
hydrogen  at  M  "  1.8  la  injected  into  a  supersonic  stream  of  air.  The 
pressure  and  temperature  of  the  air  are  sufficiently  high  so  that 
chemical  reaction  takes  place.  The  speed  at  the  center  of  the  hydrogen 
Jet  decreases  because  of  the  pressure  gradient,  however,  the  pressure 
gradient  is  small  and  therefore  the  stream  remains  supersonic.  In  the 
region  of  combustion  the  temperature  Increases  and  the  Mach  number 
decreases,  here  the  flow  becomes  locally  subsonic  and  thercfoie  more 
sensitive  to  pressure  gradients.  The  pressure  gradient  continues, 
here  the  flow  reaches  locally  zero  velocity  and  reverse  flow  is  produced 
at  some  distance  from  the  axis.  Fig.  15  gives  some  velocity  profiles 
before  reverse  flow  occurs.  When  reverse  flow  is  gene  rated , two  opposite 
vortical  rings  are  produced  (shown  schematically  in  Fig.  16),  which  are 
very  unstable.  Experimental  evidence  of  such  flow  is  shown  in  Fig.  18. 
The  flame  spreads  very  rapidly  radially  as  soon  as  reverse  flow  is 
produced  by  adverse  pressure  gradients. 
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Such  type  of  phenomena  la  of  great  Interest  In  many  practical 
problems,  for  example,  In  boundary  layer  In  presence  of  ablation 
when  the  ablative  material  reacts  with  the  air.  The  flow  for  such 
phenomena  Is  difficult  to  describe  and  analyze  because  the 
adverse  pressure  gradients  can  change  strongly  even  the  qualitative 
characteristics  of  the  flows. 

b,  CONCLUSIONS 

From  these  considerations  It  appearo  evident  that  the  problem 
of  mixing  of  heterogeneous  gases  Is  far  from  being  understood 
quantitatively.  The  analytical  approaches  available  arc,  In  many 
cases,  not  jut lsfactory ;  the  experimental  data  are  Insufficient,  and 
often  not  accurately  presented.  Many  new  phenomena,  some  characteristic 
of  the  heterogeneity  of  the  two  streams  are  present  because  of  the 
production  of  adverse  pressure  gradients,  which  even  If  small,  often 
can  change  strongly  the  flow.  Therefore, It  can  be  concluded  that 
additional  work  In  this  relatively  new  field  la  required. 
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FIGURE  6:  MATCHING  ANALYSIS  BETWEEN  BOUNDARY  LAYER  AND  FLOW  WITH  PRESSURE  GRADIENT 
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INTEGRATION  OF  FLOW  EQUATIONS 
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FIGURE  10:  MACH  NUMBER  PROFILES  AT  SEVERAL  AXIAL  STATIONS  SHOWING  THE 
EFFECT  OF  VISCOSITY  ON  FLOW  REVERSAL 
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FIGURE  17.  Theoretical  centerline  itatlc  temperature  distribution 
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